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Abstract 

Assuming a Gaussian approximation for the QCD gluo dynamics, all the nonpertur- 
bative physics can be encoded into two parameters: the gluon correlation length T g 
and the gluon condensate G2. These parameters are sufficient in order to describe 
the heavy-heavy quark nonperturbative interaction. In this work we adopt the same 
framework in order to study heavy-light bound states in the non-recoil limit. Spon- 
taneous chiral symmetry breaking and a confining chiral non-invariant interaction 
emerge quite naturally. The gap equation is solved and discussed. In particular a 
relation between the light quark condensate and T g is derived. The energy spectrum 
for the bound state equation is evaluated and commented. 
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1 INTRODUCTION 



The study of heavy-light quark bound states has absorbed in the last years a 
huge amount of effort either on the experimental side, where a lot of new data 
have been produced (compare, for instance, with [1]), either on the theoretical 
side, due to the use of heavy quark symmetries [2]. Around the limit where 
the symmetries are exact (infinitely heavy quark) a systematic expansion in 
the small parameters Aq C d/^q and q; s (Mq) hold, Mq being the mass of 
the heavy quark. This expansion contains the long-distance physics of several 
observables encoded into few hadronic parameters, which in general can be 
defined in terms of operator matrix elements in heavy quark effective theory 
(HQET) [3]. In particular in HQET the heavy-light meson mass is given by 

M q Q = M Q + A + O ( — J corrections. 
\MqJ 

The parameter A is one of the nonperturbative parameters of the HQET. It 
is the sum of all the contributions coming from the terms independent on the 
heavy-quark mass Mq. A can be fixed on the data. In order to calculate it from 
QCD one needs some dynamical input. This dynamics, being the dynamics of 
the light quark and the glue, is inerehently nonperturbative. Some approaches 
resort to calculation via phenomenological potential models [4,5], sum rules 
[6] or relativistic phenomenological equations [7]. To have a well founded cal- 
culation of A is of great importance since its value affects the determination 
of many phenomenological quantities (cf e.g. [8]). 

An attempt in this direction has been done in [9]. In that work an unified 
description of the heavy-heavy and the heavy-light dynamics in the non-recoil 
limit was suggested on the basis of the so-called Stochastic Vacuum model 
of QCD [10]. This consists in a Gaussian approximation of the QCD gluody- 
namics which is assumed to be governed by the two-point (non-local) gluon 
condensate. This assumption is confirmed by lattice simulations data [11] in 
the case of quasi-static quarks. In this way a new dynamical scale explicitly 
appears: the gluon correlation length T g . In [9] the heavy mass limit (m > T" 1 ) 
was explicitly explored. Under this condition two different situations occur. If 
T g x is bigger than all the other scales of the problem, then the one-body limit 
of the heavy quark potential is recovered with the famous Eichten-Feinberg- 
Gromes spin-orbit scalar-like interaction [12]. The other case, T~ x smaller than 
all the other scales of the problem, is the typical heavy quark sum-rule situa- 
tion. As a check the well-known expression of the heavy quark condensate as a 
function of the gluon condensate was obtained. Realistic heavy-light systems 
are characterized by m < T~ x . In this case the quark propagator cannot be 
considered as a free one and its dynamics is inherentely nonperturbative. For 
m = our effective Lagrangian is chiral symmetric. Then the most appro- 
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priate approach is to solve the non-linear bound state equation as suggested 
in [9,13] or, in a different language, to determine the physical chiral broken 
vacuum of the problem. This will be done in the present work by means of the 
Bogoliubov-Valatin variational method. Since the Gaussian approximation of 
the QCD gluodynamics works very well in the heavy quark limit [14], one 
expects reasonable results for, at least, B s and D s mesons. 

More from the point of view of principles this work would like to shed some 
light on the QCD nonperturbative dynamics. It is clear that when we deal 
with light quarks, chiral symmetry breaking and confinement are the two rel- 
evant facts. But it is also true that it is difficult to get hold of both effects. 
On one side, Nambu-Jona-Lasinio models, as well as instantons, explain the 
main physics of chiral symmetry breaking but do not have the color confine- 
ment property of QCD. On the other side, the confining models that try to 
explain the light quark phenomenology are mainly of two types: 1) models 
with a scalar confining interaction [7] that explicitly breaks chiral symmetry 
but insures the spin-orbit interaction known in the heavy quark limit [12]; 2) 
models with a vector confining interaction [15,16] which are chiral symmetric 
and allow spontaneous breaking of chiral symmetry but then are in trouble 
with the spin-orbit interaction. The heavy-light interaction which we propose 
is obtained cleanly from QCD in a gauge invariant approach and under the 
assumption of the dominance of the bilocal correlator. It has the feature of 
being able to reproduce the expected confining linear potential and the ex- 
pected nonperturbative spin-orbit coupling in the heavy quark limit as well as 
of being chiral symmetric in the chiral limit [9] . More in detail, our interaction 
kernel turns out to be more complicated than a simple convolution kernel. It 
depends on two dimensional parameters, the gluon condensate and the corre- 
lation length that control the nonperturbative physics. The string tension a 
arises as a function of these two parameters in the potential limit. We stress 
that the presence of these two parameters is relevant in order to reproduce 
the flux tube structure which is the main physical fact of the heavy quark 
nonperturbative dynamics. For heavy-light systems chiral symmetry and con- 
finement appear to be strongly related as also suggested by lattice simulations 
at non-zero temperature where confinement and chiral symmetry breaking 
show up at the same temperature [17]. This is because in the proposed ap- 
proach chiral symmetry breaking effects emerge in the physical framework of 
a quark- ant iquark bound state. As a consequence chiral symmetry breaking 
modifies the binding interaction, which turns out to be non-chiral invariant 
on the physical vacuum. This is very different with respect to the traditional 
Dyson-Schwinger methods where chiral symmetry breaking is associated with 
the (gauge dependent) non-linear dynamics of a single quark [18]. But it differs 
also from the traditional phenomenological motivated Dirac equations used in 
the literature for describing heavy-light systems in the non- recoil limit [7]. 
In those works a scalar confining interaction which explicitly breaks chiral 
symmetry is introduced by hand and phenomeno logically justified. Here a chi- 
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ral non-invariant binding interaction emerges naturally on the chiral broken 
physical vacuum. The main pieces of the puzzle seem to go at the right place. 

In the present work a simplified expression for the non-local gluon condensate 
with respect to the parameterization derived from lattice simulations [19] will 
be used [20] . This is only a technical assumption which in this exploratory work 
helps to increase the transparency of the calculation, leaving the main physics 
untouched. The main point is that we are able to give an unified description of 
the heavy-heavy and of the non-recoil heavy-light bound states in terms of few 
parameters (the gluon condensate, the gluon correlation length) determined 
by the nonperturbative dynamics of QCD. 

The plan of the paper is the following. In section 2 we set up the formalism 
and discuss with some extension the Gaussian approximation of the QCD 
gluodynamics. In section 3 we derive the gap equation and evaluate the quark 
condensate for different values of the mass. In section 4 we write the heavy- 
light bound state equation in the non-recoil limit and calculate the spectrum. 
Conclusions and outlook are given in section 5. 



2 THE HAMILTONIAN 



A system made up by a quark q of mass m and an antiquark Q of mass Mq 
is described by the 4-point Green function: 

G- mv (x,u,y,v) = (0\q(y)U(y,v)Q(v)Q(u)U(u,x)q(x)\0). (1) 

The Schwinger strings U(y,x) = P exp jig J ds {y — xYA^(x + s(y — x)) j 

have been added in order to have gauge invariant initial and final bound 
states. P stands for the path ordering of the color matrices. 

Let us consider the limit in which the antiquark is infinitely heavy (Mq — > oo). 
The heavy quark behaves then like a static source propagating from u = 
(-T/2,0) to v = (T/2,0). Let us take x = (-T/2,x) and y = (T/2,y). For 
infinitely heavy Q we have then 

G- mv (x, u,y, v) = e- m Q T JVqVql exp {i J d A x C{x)} q(y)U (y , v)U (v , u) 



xU(u,x)( ] —f] {) q(x)\, (2) 



C(x) = q(x)(ip— m)q(x), 
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(0) = —J VAe lSYM O(A), S YM = Yang-Mills SU(3) action, 

where the upperscript (2) refers to the second heavy fermion line, = <9 M — 
igA^ and M is a normalization factor. A "natural" gauge choice is that one 
which eliminates the contributions coming from the straight-line Schwinger 
strings U in Eq. (2) [21]: 

A M (x o ,O)=0, x J A j (x o ,x)=0. (3) 



We then have 



G illv (x,u,y,v) = e- lM * T JvqVq (exp J d A xC{x)^q{y) (^^) 

(4) 

In the Wilson loop language, the Wilson loop VF(r) made up by the heavy 
quark trajectory, a generic light quark path z connecting x with y and the 
initial and final point Schwinger strings (see Fig. 1) simply reduces in the 
gauge (3) to the light quark contribution only: 




Moreover with this gauge choice the gauge field A^ can be expressed in terms 
of the field strength tensor F^ u : 

l i 

A (x) = J da x k F k0 (x , ax), Aj(x) = J daax k F k j(x ,ax). 

o o 

We notice that, though this gauge appears to be natural in this problem, 
its choice is completely arbitrary and motivated only by convenience. Being 
the starting Green function (1) gauge invariant, by proper handling we would 
obtain exactly the same results within any gauge. 

In order to go on we need an assumption on the gauge field average in Eq. (4). 
Up to now this is unavoidable since the light quark sector does not provide us 
with an obvious expansion parameter. 1 We assume that 



'exp^ig J d 4 x q(x)^ A^x)q(x)^ = exp ji J d A xq(x)^T a q(x){gAl(x)) 



1 For instance, in the case of heavy quarks the system is essentially non-relativistic 
and we can expand in the quark velocity. 
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I/=(T/2,») 

x = (-T/2,x) 



V 



u=(-T/2,6) v = (T/2,6) 



Fig. 1. The Wilson loop in the static limit of the heavy quark. 

~\jd A xJ d%q(xh»T a q(x)q(y)YT b q(y)(gAl(x)gA b u (y)} 
+ higher order clusters | 
~ exp {-i J d 4 x J d 4 yq(x)YT a q(x)q(y)YT b q(y)( 9 A;(x)gAl(y)}} . (5) 

The first equality is exact. It follows from interpreting the gauge field average 
as a statistical average and from performing a cluster expansion [10]. The 
second is an approximation. It corresponds to stop the expansion with the first 
non- vanishing cluster. This approximation was used in this context in [9,13]. 
Although it has been very successful in the last years either in applications 
to the heavy quark potential where it was first proposed [10], as well as in 
the study of soft high-energy scattering problems (for some recent reviews 
see [14]). Equation (5) is our key assumption. It has been proven to work 
for quasi-static quarks as the phenomenological success of the quoted papers 
seems to show and as it is confirmed by lattice data (see discussion in [11] and 
[22]). For systems with a light quark there is no a priori reason nor to accept 
it nor to discard it. At this stage of our understanding we will assume that 
approximation (5) works and see how far we can go with it. In particular in 
this work we will see how chiral symmetry breaking emerges in a heavy-light 
bound state and how it leads to a non-chiral invariant confining interaction. 
Since we are confident that our assumption works al least for heavy quark 
bound states we will assume a conservative attitude and realistically expect 
that our picture might be able to describe B s and D s mesons. Nevertheless in 
the following we will not put a lower bound on the light quark mass. 

The effective Lagrangian density we will use is therefore 



C(x) = q(x)(i $ — m)q(x) 

+ d^q(x)rT a q(x)q(yWT*q(y)< g AXx)gJ%(y)), 
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with 



1 1 

(gA;(x)gA b u (y)) = x k y l J da a n ^ J df3 ^ (gF^(x° ,ax)gF*(y , fr)) , (6) 

o o 

where n(0) = and n(i) = 1. The quantities in Eq. (6) are gauge invariant 
since, due to the gauge choice, all the field strength tensors can be thought to 
be connected by straight-line Schwinger strings. The corresponding Hamilto- 
nian is given by: 

H — J d 3 x (x)(—ia ■ V + m/3)q(x) 

-- J d 3 x J dSq{x)rT a q{x)q{y)YT b q{y){gAl{x)gA b v {y)), (7) 

where we have introduced the traditional notation (5 = 7 and a 1 = 7°7*. 

We notice that the interaction with the heavy quark source only apparently 
disappeared from H. Actually, it manifests itself via the breaking of transla- 
tional invariance in the 4-fermion term due to the gauge choice (3). Therefore 
the effective Hamiltonian given in Eq. (7) is simply another way to write the 
bound state equation already given in [9,13]. It can, and in principle it should, 
be used as it is. The whole nonperturbative physics is contained into the non- 
local gluon condensate (g 2 F^ l/ (x)F b x (0)) and a parameterization of it is known 
from lattice measurement [19]. Being interested only in the nonperturbative 
contributions, a good parameterization is (in Euclidean space) 

(9 2 F^(x)F b pX (0)) = ^ ab (8, p 5 uX - 5, x 5 up )(g 2 F\0))e-^. (8) 

In particular it manifests a long range exponential fall off in \x\ with a gluonic 
correlation length T g ~ 0.3 -j- 0.4 fm. Once the form of the non-local gluon 
condensate is plugged in Eq. (6), the Hamiltonian describing the nonpertur- 
bative dynamics of heavy-light mesons in QCD (under the crucial assumption 
(5)) is in principle completely defined. 

At this point we are left with the usual technical problems connected with the 
actual solving of the bound state equation. In the following we will bypass all 
these by not keeping the "realistic" lattice parameterization of Eq. (8), but 
by making a rough approximation on it. This will allow us to perform almost 
all the calculations analytically and to take advantage from some existing 
literature. The approximation will be very rough and we do not expect that 
it will allow us to make quantitative predictions. Nevertheless the qualitative 
features of the Hamiltonian (7) will be preserved. In particular we will see how 
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the mechanism of chiral symmetry breaking is expected to work in the bound 
state framework and how it leads to a non chiral invariant effective bound 
state interaction. Let us, first, neglect all the perturbative contributions to 
(g 2 F^ u (x)Fp X (0)}. They are irrelevant since they do not lead to confinement 
and to chiral symmetry breaking. We assume that the leading contribution 
to it is given by the electric fields only and that they are allowed to vary 
only in time. Moreover we approximate the exponential fall off in time with 
correlation length T g with an instantaneous delta-type interaction. This can 
be done since for light quarks the energy scale T" 1 is expected to be bigger 
than the other scales of the problem. More explicitly our crude approximation 
is 

(9A a »(x)gA b u (y)) ~ J^M^ x . y S (x _ yo) Tg (g 2 E\0)). (9) 



A similar assumption can be found in [20]. The effective Hamiltonian (7) 
becomes now: 



H = J d 3 xq^x)(-ia-V + mp)q(x) (10) 
-\V 3 J d 3 x J d 3 yr 2 q\x)T a q(x)qt(y)T a q(y) 

+2V 3 J d 3 x J d 3 yR 2 q\x)T a q(x)q\y)T a q(y), (11) 
^-|(^ 2 (0)> = |- 2 G 2 , G 2 = AF 2 (0)>, 

where R = x/2 + y/2 and r = x — y. Since the interaction is instantaneous 
from now on we will consider the fermion fields as a function of the spa- 
tial coordinates only. From Eq. (11) it is very clear the role played by the 
approximation (9). It has allowed us to disentangle trivially in our effective 
Hamiltonian the self interacting part (function of r) from the external source 
interacting term (function of R). In the starting Hamiltonian (7) with a "real- 
istic" lattice parameterization of the non-local gluon condensate, like Eq. (8), 
these terms might be mixed up in a very complicate way. 



3 THE GAP EQUATION 



In this section we concentrate on the chiral properties of the Hamiltonian (7). 
Following [15,16] we will use a Bogoliubov-Valatin variational method in order 
to select the chiral broken vacuum. Let us expand the quark fields on a trial 
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basis of spinors: 



?(*) = E 



d 3 k 
J2^y 



ikx 



u s (k)b s {k)+v s (k)dl(-k) 



;i2) 



where the trial spinors u s and v s satisfy the usual orthonormality relations 
4(kK(k) = vl,(k)v s (k) = 5 SS ,, ul,(k)v s (k) = 4(k)u a (k) = 0. 

A possible choice is 



«,(k) = ^ 



\Jl + sin + ^1 — sin0(/c) a ■ k 



u s (k) 



v/2 



1 + sin0(A;) — \J\ — sin a • k 



(13) 
(14) 



where and t> Q s are the usual rest-frame spinors for the free particle on the 
chiral invariant vacuum. In a moving frame the spinors on the chiral invariant 
vacuum are given by sin0(A;) = mj\Jm 2 + k 2 and cos(p(k) = kj ' \Jm 2 + k 2 . 
In particular, in the limiting case = the trial spinors reduce to the free 
massless one, while for <f> = n/2 they reduce to free static sources. Let us also 
define the projectors 



A+(k) = ^2u s (k)ul(k) = - [l + /9sin0(A;) + a • kcos0(A;) 

s ^ 

A_(k) = ^«,(k)i;J(k) = \ [l -/3sin0(fc) - a • kcos0(A;)] = 1 - A+(k). 

The operators 6, d and &t, d' are annihilation and creation operators respec- 
tively. They define the trial vacuum. Eventually after minimizing the vacuum 
energy they will define the chiral broken vacuum. Expanding the Hamiltonian 
(11) on (12) we obtain 



h=e + h; + + H 4 



s = v 



d 3 k 



Tr[(o:-k + m/3)A_(k) 



(2tt) 3 



3 r e - i(k_k,) - r r 2 



(15) 



(16) 



(2tt) 3 J (2tt) 
H r 2 =Jd 3 x: g f (x) (-ia ■ V + m(3) g(x) : 

-\ v »l d ' R I dh ' r2 'iwf^'' gt(x) ( A +( k )- A -( k )]^) ; ( 17 ) 



8 



= f / d 3 RR 2 I d 3 r /(0e* k ' r : gt( x ) [A+ ( k ) _ A _(k)] q{y) : (18) 
H, = ~\vij d 3 R J d 3 r (r 2 - 4R 2 ) : gt( x ) T a g(x)gt( y )T a g(y) : (19) 

where : : is the normal ordering operator and V is the volume of the space. 

£ is the vacuum energy. On our trial basis it is a function of the chiral angle 
0. The Bogoliubov-Valatin variational method consists in choosing which 
minimizes £: 

5£(<f>) = 0. (20) 
This is the gap equation. From Eq. (16) more explicitly we get 



2 r 2 2 

m cos <j>(k) - k sin Mk) + -K 3 <j>(k)" + T MkY + — cos Mk) sin 0(fc) 

3 L k k 



= 0. 

f21^ 



This is just the equation solved in [16]. As noticed there, it is identical to the 
Dyson-Schwinger equation we would obtain for the self-mass of the light quark 
interacting via a harmonic oscillator potential in the instantaneous ladder 
approximation. The solution 0, which is a function of k, defines the new chiral 
broken vacuum. A plot of (f>, solution of Eq. (21), for different values of the mass 
m is shown in Fig. 2. As the mass increases as much the solution approaches 
the infinitely massive limit = tt/2. Finally, the explicit calculation shows 
that 0, solution of Eq. (21), not only is a stationary point of the vacuum 
energy, but also that his vacuum is energetically favored with respect to the 
chiral symmetric one (in the massless case £(0) — £(0) < 0). 

The light quark condensate can be calculated explicitly from Eq. (12): 

3 °° 

(0\qq\0) = / dkk 2 sin<p(k). 

7T 2 J 


On the solution of Eq. (21) we get for m = (see Ref. [16]) 

(0\qq\0) = T g G 2 x 0.3722 (22) 



The result (22) looks appealing. It establishes a connection between the gluon 
condensate and the light quark condensate. The connection is possible since 
the non-local gluon condensate has introduced into the game a finite corre- 
lation length T g . Substituting to T g the lattice value 0.35 fm [19] and to the 
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Fig. 2. Chiral angles (solution of the gap equation) for different quark masses. The 
momentum k and the mass m are in units of 4y 3 /3. 



gluon condensate G 2 the value of 0.048 GeV 4 [23] we get (0\qq\0) ~ - (110 
MeV) 3 which is a rather low value. This is not surprising and is entirely due 
to our crude assumption (9). If used for static sources Eq. (9) would lead to an 
harmonic oscillator confining potential systematically below the expected lin- 
ear rising potential with slope a = 0.2 GeV 2 in the range of interest from 0.1 to 
1 fm. This suggests that, if we would use the covariant parameterization of the 
non-local gluon condensate given in Eq. (8), that we know from Ref. [10] leads 
to a phenomenological correct linear confinement between static sources, we 
would presumably enhance the light quark condensate to the usually accepted 
value of (0\qq\0) ~ - (250 MeV) 3 . 

Let us briefly discuss the term (17). Expanding on (12) we get: 



HI 



d 3 k 



[A(k) sm(p(k) + B{k) cos0(fc)] 

x &t(k)Mk) + 4(kK(k; 
[A(k) cos (j)(k) - B(k) sin 0(Jfe)] 
x [&J(k)4(-k) + d fl (-k)& a ,(k) 



(a-k) 



ss' 1 



A(k) = m + - Vq A sin 4>(k), 



B(k) = k + - V£ A cos <j>{k) - — cos <f>(k) 

k 2 



The stationarity condition (20) or (21) cancels the so-called Bogoliubov ano- 
malous term (proportional to A cos — B sin 0) which would destabilize the 
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vacuum. On the physical vacuum we then have 

d 3 k 



H 2 = E/ J^f E ^ [ 6 !( k )W + 4(k)rf s (k)] , (23) 



with i? 2 = A 2 + B 2 . Eq. (23) simply gives the light quark kinetic energy on 
the physical vacuum. 



4 THE BOUND STATE EQUATION 



The binding is given in the effective Hamiltonian (15) by the terms and 
H4. As far as the bare qQ mass is concerned we do not need to evaluate 
i?4 matrix elements. For instance, a term like (q\HAqq) would be responsible 
for the coupling of mesons to the bare qQ state [24]. Notice that, terms like 
(OlH^lqq) which would lead to coupled channels and terms contributing only 
for baryons, have been also neglected. Expanding Eq. (18) on the quark field 
(12) and taking the matrix element between a one-particle state of momentum 
p and a one-particle state of momentum q, we have 



H*(p,q) a9l = (0|6 a (p)if 2 fl &!,(q)|0) = 



-V 3 ul(p)hsm 



p + q 



a ■ p + a ■ q 

H cos 



p + q 



X 



-A (27r) 3 5 3 (p-q) 



(24) 



As expected the binding interaction would be chiral invariant (~ ct-p+a-q) for 
a massless particle on the perturbative vacuum (0 = 0). While for a infinitely 
massive particle (0 = tt/2) chiral invariance would be maximally broken. In 
our case the solution of the gap equation (21) gives rise to a binding interaction 
which contains two pieces. One is chiral invariant and the other, proportional 
to f3, breaks explicitly chiral invariance. The existence of such a term is sug- 
gested by the spin-orbit structure of the heavy quarkonium potential whose 
relativistic origin may be traced back to a scalar confining Bethe-Salpeter 
kernel [25]. In a Hamiltonian language this would just correspond to an inter- 
action proportional to (3. Indeed such a kind of interaction has been used, also 
recently, in phenomenological applications to the heavy-light spectrum [7] . As 
already argued in [9], what one obtains following the approach discussed in 
this work is an interaction which turns out to be not only proportional to f3. It 
manifests, also under the strong simplifying assumption (9), a more complicate 
structure which interpolates between a chiral invariant vector interaction and 
a scalar interaction generated both by the mechanism of chiral condensation 
and by the initial light quark mass. 
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Substituting into Eq. (24) the explicit expression for the spinors given by Eq. 
(13) and (14) and integrating over a wave function $ we get: 



^(p,q)$(q) 



3 V ° { 2^ [1 " Sin0(p)]2 " ^ [1 " Sin0(p)] ff ' (pXV) " A J ^ 

Summing up the contributions coming from the pieces if£ and °f the 
Hamiltonian, the bound state equation on the physical vacuum reads 



{^(P) + ^O 3 (^[1-Sm0(p)] 2 

+j 2 [l - sin 0(p)] S-L-A^j $(p) = A $(p), (25) 

where we have introduced the spin operator S = a/2 and the orbital angular 
momentum operator L = r x p. The eigenvalues A of the equation are the 
energy levels of the bound state in the non-recoil limit, i.e. the difference 
between the mass M q q of the considered heavy-light meson and the mass Mq 
of the corresponding heavy quark. 

In the heavy quark limit (m = Mq — > oo, (j) 71 / 2) we get 

p 2 8 



-r Vn A 



2M Q 3 ' i; 



$(p)=A$(p). 



This is the bound state equation for two static sources. As already announced, 
assumption (9) leads to a confining potential between infinitely heavy quarks 
which is a harmonic oscillator ~ V^r 2 . 

The general spectroscopic properties of Eq. (25) are simple to derive. Since 
2 S • L = Jf — S 2 — L 2 , where J/ is the total angular momentum of the light 
quark, and A = d 2 /dp 2 + (2/p)d/dp — L 2 /p 2 , the energy levels are sensitive 
only to J[ and L. In particular for each orbital angular momentum quantum 
number I we have two different levels, one with ji — £ + 1/2 (i.e. 2 S • L = £) 
and the other with ji = £ — 1/2 (i.e. 2S-L = — 1 — £). Oi course, since in this 
treatment the heavy quark symmetry is exact, there is no dependence of the 
energy levels on the heavy quark spin. In Tab. 1 we list up to n = 4 the energy 
spectrum obtained by solving numerically Eq. (25). This spectrum exihibits 
a small deviation from the pseudo-spin symmetry (i.e. the spectrum becomes 
for J excited states, almost parity independent) which measures the extent of 
chiral symmetry breaking. The value of A Uj( j ~ 0.6 GeV which we obtain with 
the values of the parameters given in section 3, appears quite reasonable. 
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Table 1 

A in 4Vq/3 units as obtained from Eq. (25) for m = 0. J and P are respectively 
the total angular momentum and the parity of the meson. 

5 COMMENTS AND CONCLUSIONS 

In this work we have studied the heavy-light non-recoil dynamics of QCD 
adopting a Gaussian approximation for the gluodynamics. All the heavy quark 
limits (m > T g ) were already studied in [9] and give the expected results. Here 
we do not have given any constraint on the light quark mass. Chiral symmetry 
breaking and a chiral non-invariant binding interaction emerge quite naturally 
in our approach and a link is established between chiral symmetry breaking 
properties and confining interaction. In particular with Eq. (22) we establish 
a relation between the order parameter of chiral symmetry (the quark con- 
densate (0|gg|0)) and that one which in our framework describes confinement 
(the gluon correlation length T g ). A similar relation can be found in [13]. 

The actual calculations were performed under the rough approximation (9). 
This is unrealistic since it gives in the heavy quark limit a confining potential 
which is not linear . Moreover all magnetic contributions were not considered. 
This is expected to affect the levels listed in Tab. 1. Nevertheless, as discussed 
in section 2, we expect that the main features presented in this work (chiral 
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symmetry breaking, chiral non-invariant binding interaction, relation between 
the order parameters) will still hold also by using a more realistic parame- 
terization of the bilocal gluon condensate, for instance that one suggested by 
lattice simulations. Indeed, the encouraging results obtained in such a simple 
framework, not only support but also make urgent the extension in this sense 
of the present analysis. Work is in progress in this direction. 
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